Zero-Temperature Complex Replica Zeros of the +/ Ising Spin Glass 
on Mean-Field Systems and Beyond 



Tomoyuki Obuchi 3 , Yoshiyuki Kabashima b , Hidetoshi Nishimori , Masayuki Ohzeki c 

"Department of Earth and Space Science, Faculty of Science, Osaka University, Toyonaka 560-0043, Japan 
b Department of Computational Intelligence and Systems Science, Tokyo Institute of Technology, Yokohama 226-8502, Japan 
Department of Physics, Tokyo Institute of Technology, Tokyo 152-8551, Japan 



o ■ 



Abstract 



5— i Zeros of the moment of the partition function [Z"]j with respect to complex n are investigated in the zero temperature limit f3 — > oo, 
— * keeping y = fin « 0(1). We numerically investigate the zeros of the +/ Ising spin glass models on several Cayley trees and 
"^h hierarchical lattices and compare those results. In both lattices, the calculations are carried out with feasible computational costs 
by using recursion relations originated from the structures of those lattices. The results for Cayley trees show that a sequence of 
t— I the zeros approaches the real axis of y implying that a certain type of analyticity breaking actually occurs, although it is irrelevant 
for any known replica symmetry breaking. The result of hierarchical lattices also shows the presence of analyticity breaking, even 
£H In the two dimensional case in which there is no finite-temperature spin-glass transition, which implies the existence of the zero- 
ed temperature phase transition in the system. A notable tendency of hierarchical lattices is that the zeros spread in a wide region 
of the complex y plane in comparison with the case of Cayley trees, which may reflect the difference between the mean-field and 
• ^ finite-dimensional systems. 
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1. Introduction 

Disordered systems are one of the challenging problems in 
statistical physics. Especially, spin glasses have been investi- 
gated for a long time as an ideal and nontrivial problem treating 
disorder. One of the most important approaches in the spin- 
glass theory is the replica method. This method has provided 
both profound concepts and useful calculation techniques in the 
theory, which promoted the expansion of the spin-glass theory 
to other disciplines after successful construction of a descrip- 
tion of spin glasses HUSH]. 

A characteristic property of disordered systems is its sam- 
ple fluctuation of the thermodynamic quantities S 01 ■ In 
the framework of the replica method, the fluctuation, which is 
reflected in higher-order cumulants of the free energy, is essen- 
tially utilized to calculate the typical free energy. This is ac- 
tually implemented by an assessment scheme of the cumulant 
generating function <p(n) defined as follows; 

0(n)= Urn ^log[Z"] y , (1) 

where the brackets [• ■ ]y denote the average over the quenched 
disorder in the system. The typical free energy / = 
- lim^oo[logZ]y /(J3N) is obtained from <p(ri) as - {if = 
lim„^o(d(p(n))/(dn) , and any high-order cumulants can be sim- 
ilarly derived from higher-order derivatives of tp(n) with respect 
to n. 
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In the usual replica framework, to obtain the full functional 
form of 0(n), we use the analytical continuation from n e N to 
n e K (or C). This is because the exact assessment of [Z"]j for 
neKis generally infeasible except for a few solvable models. 
However, this procedure of the replica method causes a prob- 
lem: Some analyticity breaking can generally occur in <f>(n) due 
to the limit N — » oo, which is essentially incompatible with the 
analytic continuation used in the replica method. This means 
that the expression analytically continued from n e N to n e K. 
will lead to an incorrect solution for the limit n — » if the break- 
ing of analyticity occurs in the region < n < 1 . To recover 
the correct solution of <p(n), in such cases, we need to know the 
properties of the analyticity breaking of <p(n) and to modify the 
solution according to the details. This provides a motivation to 
develop a method directly investigating the analyticity breaking 
of <p(ri) with respect to n. 

This type of transitions of <p(n) is considered to be related 
to the replica symmetry breaking (RSB) in the Parisi scheme 
Jfl 01, which is considered to be exact for a wide variety of 
spin-glass models in the limit n — » (it g ives the exact result 
for the Sherrington-Kirkpatrick model Il0l.ll in t. Actually, for a 
variation of the discrete random energy model, it is shown that 
the analyticity breaking of <p(n) actually occurs and is relevant 
to the one-step RSB (1RSB) Il2lll3ll . This fact again motivates 
us to investigate the analytic behavior of <p(n) and to examine 
the phase transitions occurring in the region n > 0. 

Under these motivations, we provide a method to investigate 
the analyticity breaking of </>(«), based on the Yang-Lee descrip- 
tion of phase transitions [ 140 . In particular, we observe the ze- 
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ros of [Z"]j with respect to complex replica number n, which 
will be referred to as "replica zeros" (RZs). 

For the discrete random energy model mentioned above, this 
strategy successfully characterizes the 1RSB transition accom- 
panied by a singularity of the large deviation rate function of 
the free energy / 11511 . On the other hand for the infinite-step of 
RSB (FRSB), a possibility that the RZs cannot characterize the 
FRSB is suggested fl6Tl . according to an argument based on the 
RZs of the +J models on some tree-like systems and an analy- 
sis of the spin-glass susceptibility. These results require more 
detailed discussions about the relation between the analyticity 
of <p(n) and the RSB. 

Another interesting problem concerning the RZs is its appli- 
cation to finite-dimensional systems. It is still a subject of con- 
siderable discussion whether the RSB occurs or not in finite- 
dimensional systems. The RZs formulation can be a help to 
examine this problem. 

For a concrete progress along the direction, we here treat hi- 
erarchical lattices B171 Il8lll9ll . Although it is known that the 
RSB is absent for spin-glasses on hierarchical lattices 12^, 21 ]> 
the dimensionality of these lattices can be tuned by changing 
a parameter and also the renormalization group analysis gives 
the exact partition function. These useful properties can make 
the hierarchical lattices be a productive first step to examine the 
finite-dimensional effects on the RZs of spin glasses. 

This paper is organized as follows. In the next section, we 
briefly summarize the formulation and the results for Cayley 
trees in 11611 . In section 3, we provide a formulation and the 
RZs plots for the hierarchical lattices, and compare the results 
with those of Cayley trees. The last section is devoted to the 
conclusion. 



2. Formulation and results for Cayley trees 

2.1. Basic formulation 

Our main objective is to solve the following equation with 



respect to n; 
[Z"h = o. 



(2) 



This transcendental equation is, however, hard to solve for gen- 
eral temperatures. Then we here restrict ourselves to the zero 
temperature limit B — » oo involving n —* and Bn — » y « (9(1). 
In this limit, the relevant contribution to the partition function 
only comes from the ground state, and the RZs equation (0 
becomes 



o, 



(3) 



where E g is the ground state energy. Besides, we focus on the 
+J models whose Hamiltonian and distribution of interactions 
without ferromagnetic bias are given by 



<H 



P(Jn) 



1 / 



y \ Jus $ j' 

+ 1) + 6( J U - 1)) . 



(4) 



(5) 



This limitation restricts the energy of the system to an integer 
value, which means that eq. (f3) becomes a polynomial equa- 
tion of x = e y , and the computational cost to calculate the RZs 
equation becomes significantly reduced. 

2.2. Results for Cayley trees 

For Cayley trees, we can efficiently calculate the moment 
[Z"]jr by combining the Bethe-Peierls approach and the replica 
method lfl6Tl . Using this scheme, the RZs equation can be 
constructed and solved in a polynomial time with respect to 
the number of spins N under appropriate boundary conditions. 
However, for Cayley trees, the number of spins and the degree 
of the polynomial of [Z"]j exponentially increases as the char- 
acteristic length of the tree L (distance between the central and 
boundary spins) grows. As a result, it is infeasible to solve the 
RZs equation [Z"]j = for large L. The resultant plots for the 
k = 2-body interacting Cayley tree with the coordination num- 
ber c — 3 and the 3-body interacting Cayley tree with c = 3 
are shown in fig. Q] We can find two characteristic behavior 
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Figure 1: RZs plots for (k, c) = (2,3) and (3,3) Cayley trees. All the zeros 
lie on a line Im(v) = n/2 for the (2,3) case (left panel) but a sequence of 
zeros approach the real axis as L increases for the (3, 3) case (right panel). The 
arrow indicates the collision point expected from the analysis of the RS order 
parameter in the L — » oo limit. 

in these plots. One is for the (k, c) = (2, 3) case (left panel) in 
which all the zeros lie on a line Im(y) - n/2 and never reach the 
real axis. That is in contrast to the other one for the (3, 3) case 
(right panel) in which a sequence of zeros approaches the real 
axis as L grows. These facts imply that analyticity breaking of 
<p(n) is absent for the (2, 3) case but present for (3, 3). 

An important point is whether the analyticity breaking in the 
(3,3) case is related to the RSB or not. According to earlier 
studies, some RSB occurs in the regular random graphs, which 
are known to share many properties with Cayley trees, of the 



same parameters (k, c) = (2,3) and (3,3) 022L 12311 . This fact, 
combined with the apparent absence of the analyticity breaking 
in the left panel in fig. [T] implies that the RZs of Cayley trees 
do not reflect any RSB. Hence, the analyticity breaking of </>(ri) 
in the case (3, 3) should be considered to be a phase transition 
keeping the replica symmetry (RS). To see this, we observe the 
asymptotic behavior of the order parameter, which is given by 
the probability that the cavity field takes zero at a distance L 
from the boundary pjj$, and calculate lini£_»oo />l;0 - P* [16]. 
For the (2, 3) case, p* becomes a simple analytic function of x = 
e y . On the other hand, for the (3, 3) case, p* shows non-analytic 



Figure 2: Asymptotic behavior of the order parameter of (k, c) = (3, 3) Cayley 
tree. The solid line denotes the L — > oo solution. A finite jump of the order 
parameter p* = lim/,-^ plq occurs at x as 2.35 op 0.85. 



behavior and the result is given in fig. [2] This figure shows a 
finite jump of the order parameter at x « 2.35 « y ~ 0.85. 
This singular point is indicated by an arrow on the right panel 
of fig.Q] The approaching point of the RZs seems to agree with 
this singular point in observation by eye, which implies that the 
singularity indicated by the RZs corresponds to the singularity 
of the RS order parameter and does not reflect any RSB. 

3. Formulation and results for hierarchical lattices 

3.1. Formulation for hierarchical lattices 

In this section, we treat the +J models without ferromagnetic 
bias on hierarchical lattices. A hierarchical lattice is consisted 
from unit cells Jl7l[l8l[l9ll . The structure of the cell determines 
the dimension of the system. Here we treat a simple cell con- 
sisting of two edge spins and q inside spins (fig.EJ. Each pair of 
edge and inside spins is connected by the interaction. The con- 
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Figure 3: A picture of a unit cell of q = 3. White circles denote edge spins and 
black ones represent inside spins. 



struction of a hierarchical lattice is performed by changing a 
bond to a unit cell, and the calculation of the partition function 
is conducted by the inverse procedure. The unit cell is renor- 
malized to a bond between two edge spins by tracing out the 
inside spins. This yields the following equations expressed by 
the renormalization relations R j and Rf\ 



J't r = Rj{{Jij\), 



(6) 
(7) 



where the set {/,/} is for the bonds in the unit cell and ,, 

J 1 j 

denotes the renormalized bond between sites i and / of the 
renormalized system, and F' v ., is the free energy of the system 



consisting of all the nested spins between i' and /. Using the 
notations in fig. [3] the explicit forms for R j and Rf are 



RjiWij}) = Yjq tanrr ' (tanWs. tanhjS/y,) = J h 

i=l P i=l 

1 i 4coshtez«. 1 / i ) 
RfUJu}) = — log L - L -^-- 



(8) 



(9) 



Using these relations, we can efficiently calculate the free en- 
ergy of general Ising systems on the hierarchical lattices. 

In random systems on the hierarchical lattices, the probabil- 
ity distribution P(Jij) characterizes the behavior of the systems 
|24j,|25[]. In the current problem, however, we need the joint 
probability distribution of the bond and free energy P(Fjj, Jij). 
The renormalized distribution P'(F?f,Ji>f) is calculated from 
the original distribution P(Fjj, Jjj) as 



P-(F r ,.Jr,)- [ \dJ,,<IF n PiF n .J n ) 

J V ij 

xS (f ff - RjiUij})) 6 F' tf - F t j - R F ({F U }) 



(10) 



Again we take the zero-temperature limit /3 — > oo, n — > keep- 
ing fin — » y. In this limit, the free energy Fy becomes the 
ground state energy E,j and only takes an integer value, which is 
also the case for the bond This enables us to exactly perform 
the renormalization ( ITOb without numerical error. Once we get 
the distribution P{Ejj, J,j), the RZs equation can be constructed 
by using the distribution of the energy 
as 



(H) 



These equations ( ITOb and ( fTTT i. which enables us to exactly as- 
sess the RZs, constitutes the main result of this section. 

3.2. Results 

The procedures mentioned above, however, involve some dif- 
ficulties. For hierarchical lattices, a characteristic length L can 
be naturally defined as the number of hierarchy of the nested 
unit cells. As L grows, the number of spins increases exponen- 
tially, which leads to the rapid growth of the support of P{Ejj) 
and the exponential increase of the degree of polynomial of the 
RZs equation. These cause difficulties in evaluating the convo- 
lution ( TTOb and solving the RZs equation (fTTT i. which makes it 
hard to treat large size systems. These restrict the values of q 
and L to moderate values. For instance, in this paper, we inves- 
tigate the cases q = 2 and 3 (the corresponding dimensions are 
d — 2 and 2.58, respectively) in the ranges L < 5 for q = 2 and 
L < 4 for q = 3. The resultant plots of RZs are given in fig. [4] 

In this figure, we can find that some sequences of the zeros 
approach the real axis of y as L grows in both q — 2 and 3 cases. 
For the two dimensional case (q = 2), those would be related 
to the zero-temperature spin-glass transitions, since the absence 
of the finite-temperature phase transition are clarified by several 




Figure 4: The RZs of hierarchical lattices. The left panel is for q = 2 and 
the right one is for q = 3, In both panels, some zeros are spreading in a wide 
region, which is in contrast to the case of Cayley trees. Also, sequences of 
zeros approach the real axis in both panels, which implies that a certain type of 
analyticity breaking occur in both cases. 



current formulation becomes a useful analysis leading to further 
understanding of finite-dimensional spin glasses. 
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researches 025U26I1 . To make this point clearer, we should treat 
the limit L — » oo as the case of Cayley trees, and now analyses 
along this direction are in progress. 

Another interesting implication from fig.[4]is that the RZs are 
spreading in a wide region of the complex y plane. This is in 
contrast to the case of Cayley trees, which implies that the ana- 
lyticity breaking of <p(ri) in the hierarchical lattices might have 
different properties from those of Cayley trees. Since in gen- 
eral the continuous zeros distribution is related to the contin- 
uous singularities of the system 02711 . the RZs observed in fig. 
[4] may be related to extraordinary phase transitions, although 
they would be different from the RSB transitions because the 
absence of the RSB for the present models is strongly suggested 
I20II21I1 . This point also needs more detailed analysis and a new 
approach to overcome the computational difficulty in assessing 
the RZs is desired. 



4. Conclusion 

In this paper, we have developed a formulation utilizing the 
zeros of the nth moment of the partition function, to directly ex- 
amine the analyticity breaking of the cumulant generating func- 
tion 0(n) appearing in the replica method, and applied it to the 
+J models on Cayley trees and hierarchical lattices in the zero 
temperature limit ft — » oo, n — > keeping fin — > y « 0(1). The 
results imply the presence of analyticity breaking of the cumu- 
lant generating function <p(ri) for both lattices, even in the two 
dimensional case of the hierarchical lattices in which there is no 
finite-temperature spin-glass transition. Referring to other ana- 
lytical properties of these systems IKj 2(1 21], we can reason- 
ably conclude that those singularities revealed by the RZs are 
irrelevant to any known RSB. Direct applications of the current 
scheme to systems exhibiting the RSB should be investigated in 
the future. 

In comparison with Cayley trees, the RZs of the hierarchical 
lattices tend to widely spread in the complex y plane, which 
may be due to the finite-dimensional effect of the lattices. This 
point should also be clarified by further studies on the finite- 
dimensional spin glasses. At present, analytical approaches to 
finite-dimensional spin glasses are still rare. We hope that our 
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